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ABSTRACT 


An expert system for the thermal design of single 
phase flow shell and tube heat exchangers (STHE) has been 
developed. The system is highly interactive which gathers 
information about the problem by asking the user necessary 
questions, and giving him help and explanation if required. 

This expert system is based on logic programming 
where knowledge is seperate .from control. The inference 
algorithm used in the expert system shell is SLD-resolution. 
The query tree is developed by this resolution and depends 
upon user response. 

The design process of STHE has been studied and 
has been transformed into a knowledge base comprising of 
rules, facts, questions and explanations. The process of 
encoding knowledge into rules is also discussed. The modi- 
fied Bell-Delaware method is used for the shell side analy- 
sis and Kern's method for the tube side analysis. The shell 
side analysis considers the non-ideal cross flow as it occurs 
in an actual STHE and suitable correction factors are calcu- 
lated to take into account the deviation from the ideal cross 


flow. 



(vi) 


NOMENCLATURE 


A o 

* 

A 

^ A o^req 

A tot 

A tp 

b 

B c 

C 1 

C br' C bp 


D 


ctl 


D otl 

D s 

D t 

D w 


Constant 

2 

Total heat transfer area, equation (3.1b), m 
Constant, equation (3.12), mm"*" 1 

Actual area required for heat transfer, equation 
(3.90), m 2 

Total flow area available on the tube side, 
equation (3.76), mm 

Flow area per pass required to maintain a velocity 
v^ms"" 1 on the tube side , equation (3.77), mm 1 
Constant 

Segmental baffle cut as a percent of D g , equation 
(3.22) 

Tube field layout constant, equation (3.12) 
Constants 

Specific heat of shell and tube side fluids, 
equation (3.1a) J kg 1 K~ 1 

Tube bundle pitch circle diameter, equation 
(3.17) mm 

Outside diameter of the tube bundle, equation 
(3.37), mm 

Shell diameter, equation (3.11), mm 
Tube diameter, mm 

Equivalent hydraulic diameter, equation (3.32)mm 



\v-tx/ 


e 

F 




sbp 


w 


’ideal 


h s' h t 


J, 


J. 




Jn 


J. 


k ,k . 

S X 


Constant 

Log mean temperature difference correction 
factor, equation (3.5) 

Fraction number of tubes in pure cross flow 
between baffle tips, equation (3.28) 

Ideal tube bank friction factor, equation (3.67) 
Ratio of bundle bypass area to overall cross 
flow area S m , equation (3.39) 

Fraction number of tubes in the baffle window, 
equation (3.27) 

Ideal tube bank heat transfer coefficient, 
equation (3.60), Wx K 

Heat transfer coefficients for shell side and tube 

— '2 

side fluids, equations (3.66) and (3.87), Vtfn K 
Bundle bypass correction factor for heat transfer, 
(3.52) 

Segmental baffle window correction factor for 
heat transfer, equation (3.47) 

Ideal tube bank heat transfer factor, equations 
(3.57), (3.58), (3.59) 

Baffle leakage correction factor for heat transfer, 
equation (3.50) 

Heat transfer correction factor for adverse 
temperature gradient in laminar flow, equations, 
(3.54) and (3.56) 

Thermal conductivity of the shell and tube side 
fluid, 



(viii ) 


L bc 


l pi 

L 

s 

L sb 

L ta 

L tb 

L to 

L , 
tp 

Ltp eff 

L ts 

L wp 

V*t 

KA 

N b 


Tube wall thermal conductivity, Whf^K."" 1 

Tube bundle to shell diametral clearance, equations 
(3.13) through (3. 16), mm 

Central baffle spacing, equation (3.21), mm 
Flow direction distance between two tubes, 
equation (3. 33), mm 

Half of tube lane partition, equation (3. 32), mm 
Length of the shell, equation (3.11), mm 
Diametral shell to baffle clearance, equation 
(3.40) mm 

Effective tube length for heat transfer, equation 
(3.11), mm 

Diametral tube to baffle hole clearance, equation 
(3.43), mm 

Nominal tube length, equation (3-18), mm 
Tube pitch, equation (3. 12a), mm 
Effective tube pitch, equation (3. 25), mm 

Tube sheet thickness, equation (3. 19), mm 
Effective distance of cross flow penetration, 

(3.34) , mm 

Mass velocity of the shell side and tube side 
fluids, equations (3.44) and (3.81), Kgrrf^s” 1 
Mass flow rates of shell side and tube side fluid, 
kg s'" 1 

Total number of baffles, equation (3.20) 



(ix ) 


N c 

N ss 

N t 

N tcc 

^tcw 

N tp 

N tw 

p 

p 


Q 

o 

R 

R b 


Re s’ Re t 

Rr- , Rr 

V f t 


Number of tube rows crossed in the heat exchanger, 
equation (3.55) 

Number of sealing strip pairs, equation (3.52) 

Total number of tubes, equation (3,79) 

Number of effective tube rows crossed between 
baffle tips, equation (3.33) 

Effective tube rows crossed in one baffle 
spacing, equation (3.35) 

Number of tube side passes, equation (3.78) 

Number of tubes in the segmental baffle window, 
equation (3.30) 

Constant, equation (3.51) 

Thermal effectivess, equation (3.7) 

Prandtl number of shell side and tube side fluids, 
equations (3.46) and (3.83) 

Heat duty of the heat exchanger, equation (3.1a),W 
Heat capacity ratio, equation (3.6) 

Bundle bypass correction factor for pressure drop, 
equation (3.53) 

Reynolds number of shell side and tube side fluids, 
equation (3.45) 

Fouling factors for heat transfer, equation (3.2), 

2 -1 
rn KW x 


r^,r Q Inside and outside radius of the tubes, equation 

(3.2), mm 



(x) 


Rj Baffle leakage correction factor for pressure 

drop, equation (3.51) 
r lm Constant, equation ( 3.48) 

r ,r Constants, equation (3.49 ) and (3.52) 

R End zone correction factor for pressure drop, 

o 

equation (3.56a) 

S^ Bundle to shell bypass area within one baffle 

p 

spacing, equation (3.36), mm 

Cross flow area at shell centreline within one 
m 

p 

baffle spacing, equation (3. 25), mm 

Shell to baffle leakage area, equation (3.41),mm^ 

S Tube to baffle hole leakage area per baffle, 

p 

equation (3.42), mm 

S w Net cross flow area through one baffle window, 

p 

equation (3.31), mm 

S „ Gross window flow area without tubes in one 

wg 

p 

window, equation (3.26), mm 

S wt Segmental baffle window flow area occupied by 

p 

tubes, equation (3.29) , mm 

T Temperature, equation (3.63),K 

T , T Inlet and outlet temperatures of the fluid, 
c in c out 

equation (3.4), C 

T. T, Inlet and outlet temperatures of the hot fluid, 
h in’ h out 

equation (3.4),°C 

T . ,T^ Inlet and outlet temperatures of the shell side 
si so 

fluid, equation (3.5),°C 



(xi) 


JL * i , 

s av tav 
T ti’ T to 


*t 

T w 

U o 


V 


t 


Average shell side and tube side temperatures, 
equation (3.62), °C 

Inlet and outlet temperature of the tube side 
fluid, equation (3.5),°C 
Tube wall thickness, mm 

Tube wall temperature, equation (3.62), °C 
Overall heat transfer coefficient, equation 
(3.2), WtTV” 1 

Velocity of the tube side fluid inside the tubes, 
equation (3.77) ms"” 1 

GREEK LETTERS 


d 


Constant 


A p bi 
A p c 

A p s 
A p r 

A p s 

A p t 
A Ptotal 
A^,AT t 


Ideal tube bank pressure drop, equation (371) kPa 
Pressure drop in pure cross flow, equation 
(3.72) ,kPa 

Pressure drop in end zones, equation (3.51a)kPa 
Pressure drop associated with change of direction 
in the tube side passes, equation (3.89),kPa 
Total shell side pressure drop, equation 
(3.66a)kPa 

Pressure drop inside tubes, equation (3.89),kPa 
Total tube side pressure drop, equation (3.89b),kPa 

Absolute temperature difference between inlet 
and outlet temperatures for the shell and tube 
side fluids, equation (3.1b),°C 



(xii ) 


A T LM 

V^t 


X] 

'sw 

G ds 

e ctl 

e tp 
^ s’ 


AI 

HE 

LMTD 

MTD 

STHE 


Log mean temperature difference, equation 
(3,4) ,°C 

Absolute viscosity of the shell side and tube 
side fluid at their respective average temper- 
ature, Cp 

Absolute viscosity of the shell side fluid at 
the tube wall temperature, equation (3,61) 
Centriangle of baffle cut, equation (3, 23), rad 
Upper centriangle of baffle cut, equation 
(3.24), rad 

Tube layout angle, deg 

Density of the shell side and tube side fluids, 

. -3 

kg m 

ABBREVIATIONS 

Artificial Intelligence 
Heat Exchanger 

Log Mean temperature difference 
Mean temperature difference 
Shell and tube heat exchanger 



(xiii) 


LIST OF FIGURES 

FIGURE DESCRIPTION PAGE 

3.1 Nomenclature for shell and tube heat 23 

exchangers 

3.2 Pitch layout angle 30 

3.3 Bundle geometry 30 

4.1 Flow chart for approximate design 58 

of STHE 

4.2 Flow chart for the design of STHE 60 

4.3 Tree structure for thedesign of \6 3 

STHE 

5.1 Methodology for using the system 85 

5 

LIST OF TABLES 

Mandatory questions 90 

5.2 Iterative procedure 92 



1 


CHAPTER 1 
INTRODUCTION 

1.1 EXPERT SYSTEMS AND AI ( Sangal, *85) : 

It has been the aim of Artificial Intelligence 
(AI) to develop systems which behave like human beings. This 
means inventing machines which can interact like the human 
being, i.e. ask questions, gain knowledge and tell what they 
are doing. 

Expert systems are problem-solving programs 
that solve important problems which are generally solved by 
human experts. In this thesis, we will be dealing with an 
expert system for the design of Heat Exchangers. 

By looking at the scenario in which the 
human expert performs, one ca'n say the following about an 
expert system. An expert system, besides having the ability 
to solve problems must engage in a dialogue with the user to 
acquire the relevant details of the problems, be able to 
explain its problem solving process, be easily modifiable 
to take care of new discoveries, and be able to deal with 
partial information. 

Recently, several factors have motivated rese- 
arch on expert systems. These systems are designed to mani- 
pulate and explore symbolically expressed problems that are 
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difficult for human researchers to solve. A problem, suita- 
ble for building an expert system is one which has a number of 
possible solutions. As the number of solutions increase , 
it becomes difficult for a human being to discover the correct 
solutions. The ability of AI systems to deal with larger 
solution spaces is important in that it extends the type of 
problem that can be solved with the same conceptual tools. 

Most notably, expert systems promise to solve 
hard problems that require the best (most expensive) human 
expertise. The very codification of expertise insuitable 
form for an expert system can lead to new insights into the 
structure of the domain, or new ideas about how to teach it, 

A system is considered an expert system if it 
meets the following criteria: 

1. The system gives correct answers or gives useful advice. 

2. Interact with the user to acquire relevant data. 

3. Justify its solutions. 

4. Explain why it has asked a question, i.e., the line 
of reasoning being followed. 

1.2 STRUCTURE OF THE EXPERT SYSTEM: 

Building an expert system by writing a large 
program that incorporates all the domain knowledge prescri- 
ptively ( as procedure , program steps, etc ) does not work. 
Such a system lacks flexibility, can not deal with partial 



information easily and, is not easy to change. 

In the present expert system, the knowledge 
about the domain of the problem ( design of STHE) is diff- 
erent from the procedure of how to apply the knowledge. The 
knowledge is represented as rules and facts, and the ' expert 
system shell ’ is an intricate meachanism which decides how 
to select and apply the rules. This is essentially the'rule 
based expert system*. 

While developing an expert system, as mentioned 
earlier, enormous amount of knowledge in terms of facts and 
rules is required. The present work is an attempt to develop 
an expert system for the design of Heat Exchange Equipment. 
Srinivas (1985) initiated the work in this area, and formu- 
lated rules for the selection of Heat Exchange Equipment. 

It can be seen from the literature available on 
the design of HE that a very large number of decisions invo- 
lving logic have to be made in the actual design, e.g. deci- 
ding the type of the bundle, the type of pitch layout etc. 

Many times, some parameters may already be known to the 
designer and therefore, need not be computed. The final 
design is obtained, basically by trial and error until cer- 
tain criteria are satisfied. There is no unique solution 
to the design problem. Whichever ( solution ) suits the needs 
of the user is ultimately the chosen solution. All the 
necessary design information needs to be structured 
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properly to develop an expert system. 
1.3 PRESENT WORK; 



The objective of the present work is to develop 
an expert system for the thermal design of shell and tube 
heat exchangers, without phase change. 

1.4 ORGANISATION OF THE THESIS: 


Chapter 2 discusses the Expert System Shell 
which is the basis for the Expert System developed during 
the course of present work. Chapter 3 gives all the equa- 
tions and the procedure for the design of shell and Tube 
Heat Exchangers. Chapter 4 discusses the structuring of the 
knowledge base used in the expert system developed. Chapter 5 
discusses how to use the system and the results obtained for 
a given situation including the sample sessions. Limitations 
of the -system and suggs ,tions for future work are also 


recorded 
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CHAPTER 2 

AN EXPERT SYSTEM SHELL 

" ' r ' 1 ' 1 ’ ' IT "~ V '' 

2.1 REVIEW OF LOGIC PROGRAMING : 

2.1.1 Programing Languages 

The programing language used here is LISP. LISP 

has the following distinet advantages over other languages 

(Avron, 1982). 

a. Applicative Style- The data structure is a list. Instead 

of being described as a sequence in which operations are perf- 
ormed, Even a LISP environment is a large LISP program, 
successively. 

b. Programs as Data- The other unique characteristics of 
LISP is that the programs are represented the same as 
data. LISP program. Even a LISP environment is a large. 

c. Association - LISP symbols namely atoms can have one or 
several properties associated with different property 
attributes. 

2.1.2 Predicates, Formula, Rule, Pattern Matching: 

I. Pattern Matching: 

The database can be distinctly divided into 2 types: 

a. Pattern which can have zero or more occureances of the 
wild card and 

b. A data item or a fact which does not have a wile card 



6 


where a wild card is a variable which can take any 
value. Example of a Robot trying to dig a treasure' 

( Nilsson, 1979) Goal seeking example of MARK. 

MARK is a robot for finding out the buried 
treasure. We can define some patterns and facts for this 
example. 

(1) (At < position>) 

( Have shovel ) 

( Have map) 

(Treasure-site < position>) 

(Inview treasure ) 
are some examples of facts. 

'At 1 is a predicate which has one argument i.e. position, 

'Have' is a predicate which can have either ’ shovel ’ or ’map’ 
or any other thing as its argument. We say a pattern matches 
a fact if we are able to find a value for the wild card, such • 
that after substituting this value in the pattern, we are 
able to get the fact. e.g. if we have a fact as • 

(2) (Have shovel) 

and the pattern is 

(3) (Have ?what ) 

Clearly, the variable ? what matches the fact, and the 
corresponding match is ’shovel’ matching of a fact is always 
a corresponding fact; a pattern may match with a fact or 
another pattern also. 



7 


Pattern matching with more than or one variable can be 
illustrated by the following example: 

- 

(4) (Colour apple ?X) 

(5) (Colour ?Y red) 

The patterns (4) and (5) match the fact (3) 

(6) (Colour apple red) 

and the corresponding match of X and Y variables can be 
given as 

(7) ((? X Red) (?Y apple)) ALIST 

Alist is a list of lists containing a variable and its 
corresponding match from the fact. 

It may be noted that a variable can match with another 
variable also. 

e.g. the pattern (4) (Colour ?X ?Z) 

can match with another pattern (5) ( Colour ?Y red) giving 
an ALIST as - 

(8) ((?X ?Y ) ( ?Z red)) 

A variable is internally represented as a list ($ var X). 
This is done by the Read macro /? which automatically returns 
a list of $ VAR and the atom following ’Y*. We can also have 
patterns that can match with several facts, 
e.g. a pattern (9) (Colour ( fruit ?X) ?Y) 

can match with 

(10) (Colour (fruit. apple) red) 

(11) (Colour (fruit mango) yellow) etc. 
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II. Rule 

Now going back to our example of MARK, a robot 
to dig treasure. We define a Rule as an assertion of the type 

(12) B 4 A1,A2,A3,A4 AN 

Cond ......... N _> 0 

where A1,A2,A3 AN are called as atomic formulas Rule 

12 states that if the atomic formulae A1 to AN are true then 
B is true. 

III. Atomic Formula: 

An atomic formula is a predicate, applied to 
their respective arguments, of which at least one is a 
variable; In our example, since MARK is a treasure seeking 
Robot, he is under the control of a top level action died 
as * find treasure*. 

We can have different rules as 

(13) (have ?X) 4 ( Inview ?X) 

(14) (excavate ?X) 4—-—- (Treasuresite ?X) 

(15) (pick-up ?X) -4-™™---“ (inview ?X) 

(16) (read ?X) 4 - — ■ (have ?map) 

and fact like 

(17) (have map) 

Then the goal given as ( read ?X) will infer 
that the Robot has to ( read map) since ?X matches 
with map from the database (17) 

After reading the map, he will come to know the 
treasure site > B. 
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Applying rule (14) gives him a command to excavate the 
treasure site B. 

After exacavating treasuresite B, the ' treasure* will be 
inview. i.e. a new fact ( inview treasure) will be generated 
and applying rule 13, one can infer that he has treasure i.e. 
(have treasure). 

The above mentioned actions are done in the * Expert system 
shell*. 

IV. Types of Variables (Avron, 1982) : 

1. Open Variable- a variable which matches to any element 
of a list and binds itself to that matching value. 

2. Closed Variable - A variable may already have a binding, 
as above, and will match henceforth only to that 
particular value. 

3. Restricted Variable- A variable may have restrictions 
placed on it. These restrictions are procedurally 
attached to it, in some way, for e.g. a Boolean 
predicate must be TRUE for the variable to match. 

4. Segment Variables - They match to a sublist of any 
length, rather than to an element. 

2.2 PRESENT SYSTEM AT THE USER LEVEL: 

2.2.1 Types of Data Structures. 

The data structures in this expert system shell 
are of the following type ( Sangal, Forthcoming) 
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a. Facts- A predicate can have one or more facts stored 
as property list. A fact is defined as an unconditio- 
nally true assertion hence it is a Rule with a nil 
antecedent. 

The syntax for fact is; 

(((< predicate arg^ arg^ n >) NIL Identifier 1) 

((< predicate arg 2 j_ ...arg 2 Q >) NIL Identifier 2) 

♦ * 

• • 

• • 

* • 

(( predicate arg ml . . . . . arg mn >) NIL Identifier m)) 

where The arguments are just atoms and not variables. When 
a goal is given to the system, a fact is generated if 
rules that are applied succeed, and when the same goal 
is given once again, the system gets the values directly 
from the facts, thus saving the task of applying rules 
every time a same goal is given. Any inference is 
stored as a fact only if the context for the predicate 
is declared as True, or the inference has been made by 
asking the question. 

b. Rules : 

As seen earlier, a rule is an assertion of the form 
B A x A 2 A n 

i.e. B can be inferred if the atomic formulas 


are true 
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A rule for a predicate i_> stored as a property list under 
the following syntax; 

(((< predicate> <arg 11 > <arg 12 > . . . • ■> <arg ln >) 

( (An ) (A 12 ) A 1 ^ Identifier^) 

P 1 

( ( <predicate> <arg 21 > <arg 22 > <ar ^2n > ^ 

(( A 21 ) (A 22 ) ... (A 2 P 2 )) Identifier 2 ) 

• « 

• 9 

9 • 

( (<predicate> <arg ml > <arg mn >) 

( ( Anl ) **** ^ A m 3 ^ Identifier m )) 

m 

where at least one of the argument should be a variable. 

1 ——— > no 6f arguments 
I > no of rules. 

P = no of atomic formulas that have to 

be satisfied in order that the consequent of 
of the rule can be inferred. 

c. Context ( Avron, 1982) 

The huge data base of the system is divided 
into sets called contexts. The basic idea is to 
repace the global database with a tree of distinct 
data bases called contexts the contexts are arranged 
in a tree because each tree represent a differnt 
state of world. As this state of world changes, a 



12 


context naturally gives rise to ’descendant' contexts 
which differ slightly. Most of the information in a 
given context will be the same as in the parent context 
just above it, so to save space, only the differences are 
stored. A process at any given time uses one specific 
context as its database the current context. 

When a context for a predicate is declared 
true, the inference is stored as a fact for that predicate. 

d. Question; 

A predicate may have a question stored under 
the property list. Question is stored as a list and when 
the question has to be asked, a function 'ASK-USER* gets 
the question and prints it, and further waits for the respo- 
nse from the user. 

When the user types the response, depending upon the answer 
given, 

(1) If the response is ' Dontknow ' a tag *T' is set under 
the property DONTKNOW, and this question will not be asked 
next time. 

(2) If the response is 'Help', it prints out the various 
commands that are available, like what , Dontknow etc. 

(3) If the response is 'Why ' , it types out the rule that 
has issued the current question. 

(4) If the response is 'What', it types out the exapanation 
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available for that predicate. Explanation is usually a form 
of standard value table, or guidelines on what answer 
should be given . 

(5) If the response is none of the above, then it takes 

it as the answer to the question and then stores it under the 
property fact. 

e. Question-Proc ; 

If the question cannot be asked in a simple 
sentence, and require some procedure to be called we can 
use a question procedure. 

f. Action-Predicate: 

An action predicate when applied changes the 

database. 

2.3. DEFINING RULES AND FACTS: 

2.3.1 Facts: 

As stated earlier a fact is an assertion which 
is always true for e.g. we know that the approximate film 
coefficients and fouling resistances for different types 
of fluids are as follows: 


IXPe, 

h, W/m 2 K 

R f ,m 2 K/W 

Light liquid 

1750.0 

0.15X10" 3 

Medium Liquid 

1000.0 

0. 20 XI 0~" 3 

Heavy liquid being 
cooled 

300.0 

0.15X10" 3 

Heavy liquid being 
heated 

500.0 

0. 15X10~ 3 
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We can have predicate named App-film-foul which 
has the following arguments- 

( App-film- foul <type of the fluid> < value of film coefficient 
<value of Fouling Resistance>) 

Hence the above table can be stored as a fact using the 
syntax of facts as follows, 

(Defprop App-film-foul 


((( 

App-film-foul 

Light 

1750.0 

0.15X10" 3 ) 

nil 

1) 

(( 

App-film-foul 

Medium 

1000.0 

0.20X10"* 3 ) 

nil 

2) 

(( 

App-film-foul 

Heavy-hot 

300*0 

0.15X10~ 3 ) 

nil 

3) 

(( 

App-film-foul 

Heavy-cold 

500.0 

0.15X10 -3 ) 

nil 

4) 


) Facts) 

Liquid being heated is the cold-fluid, and liquid being cooled 
is the hot-fluid hence the type of fluid is abbreviated 
accordingly. 

2.3.2 Rules: 

We know that a rule is an assertion whose consequent 
is true only if the atomic formulas in the antecedent are 
true. One can now take the example deciding which fluid to 
place on the shell side and tube side. The following set of 
rules can be used to determine which fluid is to be placed 
on the tube side: 

The tube side will have 


1. The corrosive fluid 
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2. The severely fouling fluid 

3. The fluid having high mass-flow rate. 

These rules are in order of importance i.e. 
if none of the two fluids is severly corrosive, then test 
which one is severely fouling if both the rules fail, then 
go in for the third rule. The third rule always succeds 
since either one of the fluid has a higher mass-flow rate. 
Rules are written as: 

( Defprop tube-side-fluid 
((( tube-side-fluid ? name-h) 

(( fluid-name hot-fluid ?name-h) 

( corrosive ? name-h ?X) 

(= ?X YES)) 

ST 10) 

(( tube-side-fluid ? name -c ) 

(( fluid-name cold-fluid ? name-c ) 

( corrosive ? name -C ?X) 

(= ?X YES)) 

ST 11) 
etc. 

2.4 WORKING OF THE PRESENT SYSTEM 

The heart of this ' Expert System Shell* is the 
infer algorithm. The system accepts any query when the 
function goal is used, e.g. if the query is ( Area-req ?X) 
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then ( Goal • ( Area-req ?X) ) will try to satisfy the 
query. This 1 Expert System Shell’ has been developed by 
Dr. R. Sangal, IIT Kanpur. 

The infer algorithm takes query and Alist as 
its argument and does a series of operations as follows: 

1. It gets the order of operation for the predicate in 
the given query. If there is no order of specified, order 
of operation is taken as ( compute-pred Action-pred Facts 
Rules Question Question proc). 

2. It checks whether any properties are stored under the 
property attributes listed in the order- °P one by one. 

a. If the predicate has properties under compute predicate, 
then the arguments are passed on to the compute predicate 
after substituting for the variables in the arguments 

of the predicate. 

b. A similar procedure is followed in case of action-pred 

c. In case of facts, the facts are retrieved by the macro 

' GET-FACTS’. A function ’ INFER-UNTIL ' checks whether 
the unification succeeds with any of the facts if it 
does, it returns nil since it cannot have two bindings 
simultaneously if the context is True, Else the match 
is obtained from the fact and binding is returned in the 
Alist. 

d. When trying to apply a rule, the antecedent is also 
just like a list of querries. Each of the query is 
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to be satisfied and hence the query is changed to a new 
one from the Antecedent of the rule, and this whole pro- 
cedure is repeated. 

d. The question is printed and the user response is stored as the 
answer. 

The'=' predicate has two meanings in this system First, it 
takes 2 arguments and 

1. If both the arguments are free of variables after 
ultimate substitution, returns True or False depending 
upon whether the binding is true or not 

2. If any of the argument doesn’t have a binding but is 
a variable, a corresponding binding is created and 
returned in the ALIST thus returning value of the 
predicate as True. 
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CHAPTER 3 


METHODOLOGY FOR THE DESIGN OF SHELL AND TUBE 
HEAT EXCHANGERS 


3.1 INTRODUCTION: 

Shell and Tube Heat Exchanger (STHE) is a 
very commonly used heat transfer equipment. It has no 
moving parts; it simply exchanges heat between two fluids. 

A general impression about its design is that it is simple 
and straightforward. Although a STHE is not a sophisti- 
cated piece of equipment, yet a large number of considera- 
tions enter the design process. If one looks at the petro- 
chemical industry in past 25 years or so, it is seen that 
this industry has undergone a tremendous change as far as 
the use of heat exchangers in terms of their number, size 
and the type, is concerned. It has now become essential 
for the HE designer to be a specialist. A poor design or 
an omission of minor details in the design can lead to failu- 
res or may involve a heavy cost of repair or replacement of 
the HE. In this context, the development of the intended 
’Expert System' should incorporate the minute design details 
for the successful functioning of the HE. This Chapter deals 
with the complete methodology of the design of STHE without 
phase change. 
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The design of STHE involves the computations 
for both the tube and the shell side. The flow through a 
tube or a bank of tubes is well defined, but that over a 
bank of baffled tubes having window and manufacturing 
tolerances is not at all defined. Hence, the tube side 
calculations can be made relatively easily. For the shell 
side a large number of methods ( Devore, 1961); Dohonue, 
1955; Kern, 1950; Short, 1942; Taborek, 1983 and Tinker, 
1947, 1958 ) are available for the determination of pressure 
drop and the heat transfer coefficient. Owing to the exten- 
sive research work by Bell at the University of Delaware, 
Taborek (1983) has recommended the Bell, Delaware method, 
as the most suitable one which is used hero to develop the 
expert system. The method is described in detail in 
Taborek (1983). The standard description of STHE components 
etc. can also be found in this reference. 

3.2 LOGIC OF THE DESIGN METHOD; 

The basic criterion that given or designed heat 
exchanger should satisfy is that it should perform the given 
heat duty within the allowable pressure drop. 

The design process can briefly be summarised in 
the following 7 steps ( Taborek, 1983). 

3.2.1 Problem Identification; 

STHE are extensively used in the petrochemical 
process industry, pharmaceutical industry, chemical process 
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industries just to name a few • The process designer comes 
up with a thermal problem of exchanging heat between two 
fluids. He has the inlet as well as outlet conditions 
known from his process requirement. 

3.2.2 Selection (Fanaritis, 1976) : 

Selection is the process, in which the designer 
selects a particular type of heat exchanger for a given app- 
lication. Selection criteria are many, but primary criteria 
are type of fluids, operating pressures and temperatures, 
heat duty, and cost. 

There are an almost unlimited number of alter- 
natives for selecting a heat transfer equipment, but only 
one amongst them is the best for the given set of constraints. 
STHE are the most widely used type of heat exchangers. 

3.2.3 Selection of a Tentative Set of Parameters: 

The heat duty in KW, a heat exchanger is requi- 
red to perform is given by 

% = K C p ^ T s = ^t C pt^ T t ....(3.1a) 

The heat duty can also be written as 

Q 0 = U o A o ...(3.1b) 

where 

* * 

M s ,M t = mass flow rates of shell and tube-side 
fluids 

= specific heats of the shell-and the tube 
side fluids. 


c ,C x 
ps’ pt 
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A T g 4!iT^ = Absolute temperature differences 

between the inlet and the exit 
temperatures for the shell and 
tube side strems respectively. 


^ T LM 
F 


Total heat transfer area 
Log mean temperature difference 
(LMTD) 

LMTD correction factor 

Overall heat transfer coefficient 


The first step is to determine U 0 which is given by 


LI o = 




+ R f s + ( 


r r. 2r D r 

■f^) (rir- +<V ^ ^ 

W 0 1 1 


where 



r i’ r o 
R f s ,R ft 


w 


Convective film heat transfer 
coefficients for the shell-and 
tube-side fluids. 

Inner and outer radii of the tube 
Fouling factors for the Shetl- 
and tube-side fluids 
Thermal conductivity of the tube 
wall material 


Since the tube wall thickness is very small, the mean tube 
area is approximated as the average of outer and inner area 
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The film coefficients are assumed for the type 
of fluids and so are the fouling factors. 

The tentative design is prepared based on the 
assumed values of the film coefficients. 

3.2.4 Rating of the Design : 

A detailed analysis of the heat transfer and 
pressure drop for the two process fluids is known as Rating. 

The Rating program takes the values of flowrates, temperatures 
and pressures, configuration of the heat exchanger and fluid 
properties as inputs, and computes the required heat transfer 
area for the given duty and given type of heat exchanger 
along with the drop in pressure. 

3.2.5 Evaluation: 

The rating program gives the values of the 
actual film coefficients and the pressure drop, for both 
the streams. If these parameters are acceptable, the design 
is complete. Mostly, it is found that the required heat 
transfer area does not match with the available area, or 
the pressure drop exceeds the permissible value. Under 
such circumstances, the design has to be modified. 

If the design is acceptable one can proceed 
for the mechanical design and cost estimation. 

3.2.6 Modification of the Design Parameters: 

The assumed values of heat transfer coefficients 
seldom comply with the values of film coefficients after rating. 




Nomenclature of heat-exchanger components 


1 Stationary head— channel 

14 Expansion joint 

27 

T ie rods and spacers 

2 Stationary head— bonnet 

15 Floating tubesheet 

28 

Transverse baffles or support plates 

3 Stationary head flange — channel or bonnet 

16 Floating head cover 

29 

Impingement baffle 

4 Channel cover 

17 Floating head flange 

20 

Longitudinal baffle 

5 Stationary head nozzle 

18 Floating head backing device 

31 

Pass partition 

6 Stationary tubesheet 

19 Split shear ring 

32 

Vent connection 

7 Tubes 

20 Slip on backing flange 

33 

Dram connection 

8 Shell 

21 Floating head cover— external 

34 

Instrument connection 

9 Shell cover 

22 Floating tubesheet skirt 

35 

Support saddle 

10 Shell flange— stationary head end 

23 Packing box flange 

36 

Lifting lug 

11 Shell flange— rear head end 

24 Packing 

37 

Support bracket 

12 Shell nozzle 

25 Packing follower ring 

38 

Weir 

13 Shell-cover flange 

26 Lantern ring 

39 

Liquid level connection 


Construction types 



DOOOOOO O® CD 

* . » . ! . t » * 

o o © 

f t «• 

€> 

© CD © 






i 


m 


GOOD® 


\ \ \ * 
i \\ 

\ \ * / 
i vi/ 




■jl 


©oooo 




S52I 




r l \ 

I ' I ' ', 


to ©©©©to© 


© $ CD 


Jig 3.1 Hanenclature far the shell and tube heat exchanger 





24 


In this case, the heat transfer area can be decreased or 
increased as the case may be. Changing the HE dimensions 
has a remarkable effect on the film coefficients. 

The film coefficients can be increased by 
increasing the flow velocities, changing the baffle spacing, 
changing the number of tube side passes etc. In an attempt 
to increase the heat transfer coefficient, the pressure 
drop also increases. 

3.2.7 Mechanical Design and Cost Estimation: 

Once the thermal design is complete, the various 
parts of the heat exchanger are checked for stresses. The 
costing of the heat exchanger is done only after it is found 
that stresses in all the parts are within the safe limit. 

3.3 APPROXIMATE SIZING OF THE HEAT EXCHANGER: 


The basic design equation is 

A - 

o ~ U FAT, ■ 
o LM 


.... (3.3) 


where Q and U Q are computed from equations (3.1) and (3.2) 
Tlm» is defined by 


1. For a counter-flow HE, 


A ^lm _ 


An -F 

m 


out 


) - (T. 


out 


in 


(T 


h 


in 


'out 


In 


( T h ont - T c, ) 
out m 


... (3.4) 



2. For a single shell and tube pass HE nean temperature 
difference (MTD) is the same as that computed from equation 
(3.4), hence F=1.0. For HE with 1 shell pass and 2N tube passes, 
F is given by ( Taborek, 1983) 


(d) In [ 


e 

2EPTT+EIE 

2-P (1+R+e j J 


(3.5a) 


For a 2 shell pass and 2 tube pass HE, flow is true counter 1 
current hence 


F = l.o 


For a 2 shell pass and 2N tube passes HE, F is given by 
(Kern, 1950) 


e . (1-P) 

2T R-1 ) 1 (l-PR) 

F = — — — 

2/P-1-R+ (2/p ) YTT-PTT l-PR ) +e 

In [ — — 1 

2/P-1-R+ ( 2/P ) VTT"P)( i-PR ) -e 

where 


R 


Heat Capacity Ratio 


T si~ T so 

* T to“ T ti 




• . . . ( 3 . 5b ) 


.... (3.6) 


P 


Thermal Effectiveness 
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T'si'Tso = Inlet and outlet temperatures for 
the shell side fluid 

T ti ,T to = Inlet and outlet temperatures for 
the tube side fluid 

5 = Inf(l-P JT (1-PR) J Ir £ 1 

= il£[ 

p 1 R 1 


The MTD is computed as 


M 


f At 


LM 


.... (3.9) 


3. For computing U Q from equation (3.2) suitable values 

of h s » h t ,R f s ,R ft are assumed * 

4. The total area required is computed using equation (3.3) 

5. The surface area, A , can be written as 

o 

A q = A q (D s ,L s ,D t , type of pitch lay out ) .... (3.10) 

where 

D = Diameter of the shell 

V 

L g = Length of the shell 

D^. = Outside diameter of the tube 

The type of pitch layout is characterized by an angle ©. 

*CP 

shown in Fig. 3.2. The ratio L /D (Called Aspect Ratio) can 

o o 

lie anywhere between 3.0 to 15.0. Larger value gives a smaller 
size HE hence economical. Generally an aspect ratio of 8.0 
for the shell is remommended . The relation between the heat 
transfer area, A , and other parameters of equation (3.10) is 
given by 
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where 



(I0" 6 )x A*[L ta (D ctl ) 2 ] m 2 ...(3.11) 

Length of the tube between two tube-sheets 
( taken equal to the shell length, L g , for 
approximate design) 


jj c -|.]_= Tube bundle diameter (taken equal to the shell 
inside diameter, D , for approximate design) 

A is the tube layout density parameter defined as 

A = 0.78te ( ^r~) 10111 ...(3.12) 

U 1 ^ L tp ; 

L^ = Tube pitch = 1.25 

= Tube field layout constant 

= 1.0 for 0^ = 45° or 90° 

= 0.866 for 0^.p= 30° 

0 = Tube layout angle (see Fig. 3.1) 

*tp 

The shell diameter ( D s = D c ti) and the len 9 th ^ L s =L ta^ can be 
computed from equation (3.11) once the aspect ratio is known. 


3.4 SHELL SIDE PARAMETERS : 

3.4.1 Shell Dimensions 

Knowing L, and D , the shell dimensions can be 

t3 S 

calculated as follows: 

1. Bundle Shell Clearance (L bb ) 

A suitable tube bundle is selected based on the user's 
requirement and the bundle shell clearance is calculated 
based upon the equations given below. 

a. Bundle type • = fixed tube sheet 
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L bb = 12.0 + 0.005 D g mm ....( 3.13) 

b. Bundle type = U tube sheet 

L bb = 12.0 + 0.005 D g mm ....(3, 13a) 

c. Bundle type = Packed Lantern ring 

L bb = 25,0 + °* 0175 D s mm ....(3.14) 

d. Bundle type = Outside packed floating head 

L bb = 25.0 + 0.0175 D s .... (3.14a) 

e. Bundle- type = Split backing ring floating head 

L bb = 25.0 + 0.175 D s mm (3.14b) 

f. Bundle-type = Pull through floating head 
inlet pressure for shell < 1000.0 kPa 

L bb = 80.0 + 0.0325 D s mm ....( 3.15) 

inlet pressure for shell j> 1000.0 kPa 

L. k = 80.0 + 0.0413 D mm ....(3.16) 

bb s 
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3.4.2 Baffle Geometry; 

1. Baffle Spacing (L bc )-A uniform baffle spacing is 
assumed initially, equal to E s ~which is used to 
compute the Number of baffles as 


N b -hj 1 * - 1 ) 

s 


(3.20) 


This is rounded off to the lower integer and the exact 
central baffle spacing is then calculated by 


'be “ N b +1 


(3.21) 


2. Segmental Baffle Cut, B , as a percent of D . 

v O 

The value depends upon the Ratio i- bc /D s and is given by 


B„ = 16.25 + 18.75 (L, /D ) 

C DCS 


..... (3.22 


Centriangle of Baffle Cut (©^ ) is the angle subtended 
at the center by the intersection of the baffle cut 
and the inner shell wall, see Fig. 3. 3 ). It is given 


2 cos" 1 [1-2 (j^§- ) J rad 


(3.23) 


Upper Centriangle of Baffle Cut (©c-^) - * is the angle 
substended at the center by the intersection of the 
baffle cut and the tube bundle diameter, (see Fig. 3. 3) 


It is given by 


, D 2B 0 

2 cos [ ( 1- (jo§:))3 rad 

C"tl 


(3.24) 



Tube layout geometry basic parameters 





90 ° Ltp L tp 



Fig 3.2 Pitch layout angle 



Pig 3.3 Bundle Geaeetry 
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3.4.3 Bundle Type ( Lord, 1970, ± 93 ® ana Kern, 1950); 

1. Fixed tube sheet HE has the simplest type of shell 
design. The tube sheets may be welded or bolted to the 
shell. Since there are no packings or gaskets involved, 
this type of HE has the maximum possible protection against 
leakage of shell side fluid to ambient. The fixed tube- 
sheet HE does not have any provision to accomodate the 
differential expansion of the shell and the tubes. 

2. Removable Bundle Exchangers; 

Shell and Tube Floating Tubesheet 

It has straight tubes secured at one fixed 
end and the other movable end is known as floating tube 
sheet. 

(i) Outside Packed Staffing Box- The shell side fluid 
is sealed by rings of packing compressed within a stuff ing 

by a follower ring. It can be used for pressures upto 
2000.0 kPa and temperatures upt 300°C. They are not used 
when even slightest amount of leakage to the ambient cannot 
be tolerated. 

(ii) Outside Packed Lantern Ring - The shell side and 
tube side fluids are each sealed by seperate rings of 
packing seperated by a lantern ring provided with weep holes 
so that leakage through either of them should be to the 
outside. Number of tube side passes is limited to one or 
two. They are used for shell side pressure < 1034. kPa 
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o 

and temperatures on shell side < 260.0 C . They have a 
direct advantage that the tubes can be cleaned very easily 

(iii) Pull Through Bundle - This type of HE has a seperate 
head bolted directly to the floating tube sheet. Both the 
assembled tubesheet and head are small enough to slide 
through the shell, and tube bundle can be removed without 
breaking the joints at the floating end . 

Though this type has an advantage of very low 
maintenance time, it accomodates the smallest number of 
tubes for the given shell diameter and the bundle bypass 
areas are larger than all other types. 

(iv) Inside Split Backing Ring- The inside split backing 
ring type overcomes the disadvantages of Pull through 
bundle. It differs from the Pull through type in the use 
of a split ring assembly at the floating tubesheet and an 
oversized cover which accomodates it. Clearances between 
the outermost tubes and the inside of the shell are same 
as that for outside peaked stuffing box type HE. 

3.4.4 Various Flow Areas for the Shell Side Stream: 

To compute, the correction factors for heat 
transfer and pressure drop , calculation of various leakage 
and flow areas is essential. 

1. Cross flow area at shell centreline within one 

baffle spacing , S m , is given by 
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where 


(3.25) 


Sm = L bc C L bb + L tp '“ f < L tp- D t^ 


Ltp Q ff = effective tube pitch 

= L^.p for 30° and 90° layout 
= 0.707 x for 45° layout 


Baffle Window Flow Areas 


Gross window flow area without tubes in the window is 


given by 


2 / ®ds 


c — (n )•*- ( a 
5 wg 4 ^ u s' K 2% 


sin ©, 


(3.26) 


Fraction of number of tubes in the baffle window. 


F , is given by 


sin 9 


(3.27) 


The fraction of number of tubes in pure cross flow 
between baffle tips is given by 


F c = 2F w 


.... (3.28) 


The segmental baffle window flow area occupied by, 


N , tubes is given by 
S wt = N t F w < I 


(3.2() 


The number of tubes in one window is given by 


N. = N. F 
tw t w 


(3.30) 
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Hence the net cross flow area through one baffle 
window S is 


w 


^v\,g ~ S wt 


(31) 


3. Equivalent Hydraulic Diameter: 

The equivalent hydraulic diameter is required 

only for the pressure drop calculations in laminar flow, i.e. 

if R < 100. It is calculated from the classical difinition 
es 


of hydraulic diameter and is given by 

4S. 


D. 


w 


w 

TtD t N tw + (it Dl 


(32) 


4. 


s ds' 

Number of Effective Tube Rows in Cross Flow: 


The determination of effective number of tubes 
in cross flow is essential for the calculation of heat 
transfer coefficient and pressure drop, and the corresponding 
correction factors. 


The number of effective tube rows crossed 
between baffle tips, N t , is ^iven by 

D_ B 


N. 


tee 


f 1-2 ( — £■ ) 1 
L L x z K 100 ;J 
PP 


(3.33) 


where 


L = Flow direction distance between 2 tubes 
PP 

= 0.866 L t for G p = 30° 

= L tp for ®p = 90 ° 

= 0.707 L t for 0 p = 45° 

The determination of the number of effective tube rows crossed 
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in a baffle window, M^ cw , depends on the flow pattern. 
Highest flow velocity exists just below the baffle tips 
and then decreases rapidly. Finally, the effective distan- 
ce of cross flow penetration can be determined by 



0.4 


[D s< 


h X 
100 ; 


D s- D ctl 


(3.34) 


This distance, L w p, is crossed twice in a baffle window, 
and hence effective number of tube rows crossed is 


N. = 
tew 


0.8 r n { B c x D s" D ctl 
L L s 100' “ 2 

PP 


] 


... (3.35) 


5. Bundle Shell Bypass Area Parameters: 


The flow area between the shell inside surface 
and bundle outside diameter is the main area where flow 
can bypass the desirable path through the tube field. The 
flow in this bypass lane can attain considerable order of 
magnitude (20-30^), thus decreasing heat transfer and press- 
ure drop. The bypass area within one baffle spacing is given 


by 

S b = 

L bc [< D s- D otl> + L pl ] 

.... (3.36) 

where 

D otl = 

Outer diameter of the tube 

bundle 


== 

D ctl +D t 

.... (3.37) 


kpl “ 

1/2 of tube lane partition, 

(mm) 


= 

D t /2 

.... (3.38) 


For calculations of the correction factors and R^, the 
ratio F^p, of bundle bypass area to the overall cross 
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flow area is required 

‘ F , = S./S m 

* # s op b m 


... (3.39) 


6. Shell to Baffle Leakage Areas? 

This area is required for the calculation of the 
correction factors for baffle leakage effects and R^. This 
area calculation involves diametral shell to baffle clearance, 
given by 

L sfc) = 3.1 +0.004 D g ...(3.40) 

The shell to baffle leakage area is given by 

L , 271-“© ■, 0 

S sb = * D s ( 2 ^) ( — 2^T £) h™ 2 ....(3.41) 

7. Tube to Baffle Hole Leakage Area (S.^): 

This area is also needed for the calculation of 
correction factors and R^. It is given by 

S tb = ? C (D t +L tb )2-D t ) - D t ] N t (1 ~ F w> (3 ‘ 42) 

where 

L^ b = Tube to baffle hole clearance 

= 0.8 mm for calculations ....(3.43) 

3.4.5 Dimensionless Parameters: 

1. Shell Side Reynolds Number - The shell side mass 

velocity is given by 

M , Q 

i?L = x (10 D ) kg/m Sec 

S 


(3.44) 
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The shell side Reynolds number is given by 

«% = D t / ltl s / t) s (3.45) 

where n = Absolute viscosity of the shell side fluid 

o 

at average shell temperature in centipoise. 

2. Shell Side Prandtl Number : 

This is given by 

Pr s - < c ps V k s> * < 10 ~ 3 > (3 ' 46) 

where k = Thermal conductivity of the shell side fluid 

o 

3.4.6 Correction Factors: 

The heat transfer coefficient for ideal cross flow 
over the tube banks ^ can be readily found out using equation 
(3.60). In an actual HE, the flow deviates from the ideal 
case. The various correction factors take into account this 
deviation. The value of each correction factor lies between 
0 and 1 for the ideal case. When all these correction 
factors are multiplied together with the ideal heat transfer 
coefficient, one gets the actual heat transfer coefficient. 
Same holds good for the pressure drop also. 

1. Segmental Baffle Window Correction Factor (J ) - This 

V 

factor considers the effects of baffle window on heat 
transfer factor which is based on pure cross flow. 

J reaches a value of 1.0 for baffle cuts around 
c 

25^ and even more than 1.0 for smaller baffle cuts. This is 
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so because of the fact that i s computed at the largest 
cross flows section, but as baffle cut decreases, flow 
velocities increase. This in turn is compensated by the 
fact that fewer tubes exist in the baffle window. 

The net effect of these is of significance. An 
approximation of the Delaware method curve is given by, 

J = 0.55 + 0.72 F (3.47) 

V V 

2. Correction Factor for Baffle Leakage for Heat Trans- 
fer,!-^, and Pressure Drop, - There is a significant 
amount of pressure difference between the two adjacent 
compartments. This pressure differential forces the fluid 
through two leakage areas- (l) shell and the baffle circu- 
mference (2) tube and the baffle tube hole. This decreases 
the effective cross flow stream and consequently h and p 
(total shell side pressure drop). These leakage streams 
can reach a considerable order of magnitude (40 yi) and hence 
are of great importance. Amongst the two leakage streams 
considered, shell to baffle stream is most detrimental to 
heat transfer as it does not exchange heat with the tubes. 

The tube to baffle leakage stream partially 
exchanges heat with the tubes hence not as disastrous for 
the heat transfer. 

The following parameters are computed to calculate 
the heat transfer and pressure drop correction factor! 
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■lm 


and 


S sb + S tb 
S sb 

s~+~57““ 


.... (3.48) 
(3.49) 


sb tb 

The correction factors for heat transfer, and for press- 

ure drop, R^, have the following characteristics: 


a. Most severe condition occurs for r =1, which correspon- 

s 

ds to the case of all leakage taking place in the shell 
to baffle area. 


b. The least severe condition corresponds to r =0, when 

s 

all the leakage is through the tube-baffle holes. 

A well designed HE should have values of not less 
than 0.6. The possible remedies for increasing are: 

a. Wider baffle spacing which will shift r^ m towards higher 
value. 

b. Increasing tube pitch or changing tube layout to 90° or 
45° will also have a similar effect. The correction 


factors 

are computed from the following 

equations : 

J i = 

0.44 (1-r )+[l~0.44 (l-r ) ] exp 

o o 

(-2.2 r. ) 
lm 

(3.50) 

h - 

exp [~ 1.33 (l+r s ) r lm P ] 

• * o *> o ( 3 • ) 

where p = 

[0.15 (l+r )+0.S] 

o 



3. Correction Factors for Bundle Bypass Effects for Heat 
Transfer, ( , and Pressure Drop (R^) • 
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The flow resistance along the gap between tube 
bundle and the shell inner wall is substantially lower than 
that of the tube field. Naturally, the stream tends to 
flow through this flow area. 


U- tube sheet and fixed tube sheet bundles have 
a small bundle shell clearance, hence the effect is smaller. 
In case of pull through bundles, this gap has to be blocked 
by sealing strips. Sealing strip are used in pairs ( on 
either sides) when exceeds 30.0 mm. 

The correction factors and R^, are given as 

follows : 


J b = exp (- 


C bh F S bp f 1 - < 2 r SS > 1/3 ^ 


(3.52) 


with a limit 


1 at r > 1/2 
ss — 


where 


and 

where 


'bh 


= 1.35 R e < 100 

s *“ 

= 1.25 .....Ite < 100 

s — 

N 


ss 


N. 


ss 

r~ 

tcc 


N gs = no of sealing strip pairs 




exp [- C bp F sfep (l-(r ss ) 1/3 )] 


(3.53) 


where 


1 = 4.5 for Re < 100 

op $ 

= 3.7 for Re > 100 

s’” 


and other parameters have the same meaning as before. 
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4. Heat Transfer Correction Factor for Adverse Temperature 
Gradient in Laminar Flow (j ), 


The Delaware data on laminar flow (R <20) exhi- 

es 

bited a large decrease in heat transfer which is eventually- 
postulated as an effect of adverse temperature gradient 
developed through the boundary layer. 

As Re,, increases, momentum change or crtial 

o 

effects begin to disturb the laminar boundary layer until 
it almost vanishes at Re = 100. 

Since the ideal tube banke curves are based on 
10 rows of tubes, the corresponding correction factor can 
be expressed as 


T ( 10 >0.18 1.51 , Q 

J r = ( TT”"' = HO 0 - 18 

c c 


(3.54) 


where 


N = no of tube rows crossed in the HE given by 


N c = ( N tcc +N tcw) ^b +1 ^ 


(3.55) 


For Re > 20 to Re = 100 
s — s 


1.51 

(N c ) 


0.18 


20- R 

+ r ®£) r 

+ v 80 Jl 


1.51 


(N c ) 


0.18 


-1 ] ...(3.56) 


Minimum value of J r = 0.4 
For Ro >100, J = 1 

S -T 


5. Correction Factor for Pressure Drop in End Zones: 

For equal baffle spacing, this accounts for the 
pressure drop at the inlet and outlet baffles and is 
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R g = 1.3 for U- tube bundles 

= 2.0 for other cases ...(3.56a) 

3.4.7 Actual Heat Transfer Coefficient on the Shell Side: 


This method is based on j. and f- factors from 

x 

the data on ideal tube bank. The ideal tube bank factor 
is written as 


J'i = jj_ (Re 

For a fixed value of 
for j ^ are 




D. 


^tp/^t 


1.25, 


the fitted equations 


©^p = 30° staggered 

j. = 1Q [0. 0533 x 2 -0. 7965 x + 0.1999] ....(3.57) 

©tp = 45° staggered, 

1Q [0.0576 x 2 - 0.8282 x + 0.3 ] (3.58) 

“*i 

©^.p = 90° inline, 

1Q [0.05386 x 2 -0. 75432 x] .... (3.59) 

^i 

where x= log-^Q ^ G s ^ 

The ideal tube bank heat transfer coefficient then becomes. 


h 


-2/3 


s ideal~ J i C ps A s ^s ^ 


. . . . (3.60) 


where is determined from one of the equations given qbove. 

((J) ) r = viscosity correction factor that accounts for 
the viscosity gradient at the tube wall versus the viscosity 
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at the average bulk fluid temperature, rj . The term ((j) ) r 

O v> 

is computed as follows: 

1. For liquids - 


«t s ) : 


n 0.14 

= (-§— ) 
v ri J 
'sw 


(3,61) 


where rj is the shell fluid viscosity determined at the 

SW 4 

J 'vbe wall temperature, T 


w 


and 


(j) s > 1.0 for shell fluid heated 
(]) s < 1.0 for shell fluid cooled. 


In order to determine n , it is essential to determine T 

SW w 

which is estimated as follows using the approximate values 


of h^ and h g 


T s “ T t 

T - T 4. ( av av >1 

w t av i+Chf/bs) 


(3.62) 


where T T, denote the average shell and average tube 

S A ,,y t, 


av* ’av 

temperatures, both of them being the arithmetic means of 
inlet and outlet temperatures respectively. It may be noted 
from equation (3.62) that the tube wall temperature approaches 
the temperature of the fluid with higher ’h’. 


If viscosity of the fluid is known at two temper- 
atures other than T and T a relation of the from 

w s y av 9 

r) = aT^ ....(3.63) 

can be used for extrapolating the viscosity within reasonable 
temperature limits. The temperature, T in equation (3.63) 
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is in degree Kelvin. Once, r) and ri are known, (d) ) r 

7 s 's,w 7 VY s 

can easily be computed. 


2. For gases, the viscosity is a weak function of temperature 
and the correction factor is formulated as 


For gas being cooled : (4> ) r = 1.0 .... (3.64) 

I +273.15 
s 

For gas being heated : (^ s ) r = ( x~HtTT5” ^ °* 25 •••• ( 3 ‘ 65 ) 

w 1 " * 

for a gas being heated, T is always higher than T and 

w ' 3 s , • av 

hence ((f) ) r < 1.0 

The actual heat transfer coefficient is calculated as 


h = h (J .J,.J ,.J ) .... (3.66) 

s s ideal c i o r 


3.4.8 Shell Side Pressure Drop ( K P ); 

s 


The shell side pressure drop is composed of three 

distinct components , A p„ for pure cross flow, A P for all 

c w 

baffle windows and Ap c for end zones. Total pressure drop 
Ap s =Ap c +Ap w + Ap c (3.66a) 

1. Ideal Tube Bank Pressure Drop : 


The ieal tube bank friction factor is given by 


= (io 3 ) ($ ) : 

2 (m s ) 2 s 


(3.67 


= f (Be ) 
s 

The fitted equations for f^ are 
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,[0.01745 x -1.1516 x + 2.162] 


... ( 3 . 68 ) 


45°, 


[0.01535 x-1.362 x + 1.821] ... (3.69) 


10 [0. 099089 x 2 - 0.9826 x + 1.2852] .... (3.70) 


where x = log^g(Re ) 


The ideal tube bank pressure drop is readily computed as 


o (m ) 2 

AP bi = 2xl0-y fl x N tcc x y- (^) s r r 

s S 

3 

where, ^ = density of shell side fluid, kg/m . 


(3.71) 


2. Pressure Drop in Pure Cross Flow pc is given by 


APC =Ap bi (N b -1) (R b ) (R A ) 


(3.72) 


3. Pressure Drop in All the Baffle Windows Crossed ( P ) i 

Mass velocity based on geometric mean of the cross 
flow area S m is given by 


- x 10^ kg/m 2 Sec 


.... (3.73) 


The pressure drop in all windows crossed, P w is given by 


for Re s > 100, 
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„ , 26(l V''’ S r N tcw . L bc 

b 1 L tp- D t (D„) 




tew be ro.m”3 w I’jp vp~ 

l“7d: + ,tt 2 +[2xl ° Tr J)R e; ,kP 


Pressure Drop in the End Zones is given by 


•’Ss "S 
. . . . (3.75) 


A P, 


• P bl < 1+ > R b R s 


(3.75a) 


3.5 Tube Side Parameters (Kern, 1950); 

3.5.1 Tube Side Passes; 

The number of tube side passes is determined on 
the basis of the flow velocity in the tubes. The flow velocity 
should be between 1.0 to 3.0 ms A The lower limit is to 
overcome fouling and the upper one to reduce corrosion and 
vibration. The total flow area available is given by 


(Dt-tf) x Nt 


(3.76; 


The area per pass required to maintain a velocity 


v_|_ through the tubes is given by 
& t xl0 6 


St v t 


(3.77) 


where , 


= thickness of the tubes, mm 
= density of tube side fluid, kg/m" 


The number of tube side passes are 
N. _ = Integer ( ) 


(3.78) 
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3.5.2 Total Number of Tubes is given by: 



N t 

= 0.78 D ctl 2 

C h L tp) 2 


where 

5 

= tube layout constant 

for 

®p 

= 30°, C 1 = 0.866 


for 

®p 

= 45° and 90°, C 1 = 

O 

« 

t — i 

3.5.3 

Mass 

Velocity in Tubes is 

given 


*t 

= M t xl0 6 

(WV 



(3.79) 


(3.81) 


3.5.4 Tube Side Reynolds Number (Re.^ 


Re t “ 


( D t — t t )A t 


(3.82) 


where r|^ 


viscosity of the tube side fluid, at average 
tube temperature, cp, 


3.5.5 Tube Side Prandtl Number (P ) 

r t 


It is given by 


Pr 


t 


lx 10 


-3 


(3.83) 


where k , 


Thermal conductivity of tube side fluid, 


3.5.6 Actual Tube Side Heat Transfer Coefficient, h^ 
Sieder and Tate's relation gives 
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h,0 A w r), 0.14 

= 1.86 ( “ rf) (r 1 ) Re < 2100 .... (3.84) 

*w z 

= 0.027 (Re t ) 0o8 (p r J 1 / 3 ((|)) 0 * 14 Re t >2100 

k 1 

(3.85) 

The heat transfer coefficient is computed from the curves, 

. _ •, n [-0. 2624 x 2 + 3.273 x - 7.412] (3.86) 

J h - xu 

where x = log 10 (Re t ) 

The actual heat transfer coefficient is then given by 

h t = j h ^ ( Pr) 1/3 ((j)) r ......(3.87) 

"t 

where (cj) ) r is to be computed for the tube side fluid in 
a manner similar to section 3.3.7. 

3.5.7 Tube Side Pressure Drop; 

1. Tube side friction factor is given by Sider and Tate's 
correlation for fluids being heated or cooled, in form of a 
graph (Kern, 1950). 

A straight line fitted to the graph is for 

Re t < 1000, 

f t = 144 ( lo(~ Log 10^ Re t)~ 0 * 3 )) .... .(3.88) 

and for Rg^> 1000, 

■ f . = 144 ( 10 [ -°' 25 (Log 10 (Re t ))-2.55] ) 
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2. The Total Pressure Diop is Composed of 2 partsjviz, 
the pressure drop inside the tubes, p^ and the pressure 
drop, p , associated with the change of direction in the 



where N. 


No of tube side passes ; 


The total tube side pressure drop is 

AW = AP * 


.. (3.89b) 


3.6 The Actual Overall Heat Transfer Coefficient U Q is computed 
from equation (3.2) and ”sing the values for h g and h^ from 
equation (3.66) and equation (3.87) respectively. Using this 
value of U , The actual area required for the heat transfer 
is obtained from 

(A) = ...... (3.90) 

' o req U Q F (LMTD; 

Area required should always be less than the area available, 
while the pressure drops on the shell side and tube side are 
within the permissible limits. A well designed HE should bo 
able to perform the heat duty by fully utili zing the allowable 
pressure drop on shell and tube sides. 
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CHAPTER 4 

RULE SETS FOR THE EXPERT SYSTEM 

4.1 FLOW CHART FOR THE DESIGN : 

% 

The overall design process of the STHE is 
quite complicated and lengthy, hence it was necessary to 
break down this process into distinct blocks and develop 
the system accordingly. The main blocks identified were, 
(Refer Fig. A .2) 

1. Approximate design 

2. Evaluation of geometric parameters 

3. Correction factors for heat transfer and pressure 
drop 

4. Actual shell side hoat transfer coefficient 

5. Total shell side pressure drop 

6. Actual tube side heat transfer coefficient 

7. Tube side pressure drop 

8. Comparison of the results and iteration 

The flow chart and the tree structure is 
presented in Fig. 4.2 and 4.4 respectively. Flow chart 
is the representation of transfer of control whereas the 
tree structure is the structure generated in trying to 
satisfy the querries to acheive the goal. 

4.1.1 Detailed Analysis of Approximate Design: 

The approximate design is a tentative set 
of heat exchanger parameters and if the design is accepted 
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No. 



after rating then this becomes the final design. This' step 
can be further broken down into ( Refer Fig. 4.1) 


(i) Compute overall heat transfer coefficient 

(ii) Compute heat rate required 

(iii) Compute the area of heat transfer required 

(iv) Design the geometry 

To compute the overall heat transfer coefficient , we need 
a tentative set of film coefficients. The film coefficients 
can be selected based on the type of the fluid. Since the 
approximate film coefficients for different type of fluids 
are known, we can store these values as facts. The overall 
heat transfer coefficient requires determination of shell 
side and tube side fluid. 


At this point we should take the decision 
regarding the type of data structure we are going to use. 

Here there can be two types i.e. for different 
properties like inlet temperature, outlet temperature, sped 
fic-hoat, density etc, for both the fluids 


(A) 

(fluid-1 

density 

<val>) 


(fluid-1 

inlet temperature 

<val>) 


(fluid-1 

outlet temperature 

<val >) 

(B) 

(density 

f luid«l 

<val>) 


(inlet temperature fluid-1 

<val>) 


(outlet- 

temperature fluid-1 

<val >) 
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Type A is centred around the fluids (fluid-1, fluid-2, etc.) 
while type B is centred around the properties (density, 
inlet temperature etc.) . Since there are a small number 
of fluids viz, 2 type A has the advantage that the ques- 
tions, explanations etc. are reduced in number. But at 
the same time, it reduces the flexibility of asking differ- 
ent questions, and also giving different explanations. The 
data structure type B does not have such drawback since each 
question and explanation can be different for each property. 
It also has another advantage that we can refer to fluid-1 
and fluid-2 as shell-side-fluids and tube-side-fluids. 

Once the data for the fluids is known, the 

overall heat transfer coefficient is computed. Computing 

heat rate is a straightforward process. The area of heat 

transfer required is a function of heat rate. Log mean 

temperature difference, overall-heat transfer coefficient 

and LMTD correction. Once the area is known, we have a 

choice as there can be any combination of shell diameter 

and shell length for the given heat transfer area . This 

problem can be solved by either keeping the aspect ratio 

fixed or keeping the shell diameter fixed. For the given 
shell diameter and length, we have a set of rules for 

designing the different components and parameters of the heat 

exchanger. The heat exchanger geometry should be modified 

to accomodate the various process constraints. This 

concludes the first step in design process. 
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2. After this brief discussion about the various steps 
involved in the approximate design, we shall see the actual 
predicates, actual rules and actual facts used in this 
system. 

a. The aim is to determine the shell diameter and shell 
length for the given heat duty. 

ST43; (Shell-Diameter ?D-S) 4-—- (Area-Zero ?A-0), 

( Area-* ?A-*), (Aspect-Ratio ?L-By-D), 

( = ?D-S (** (•* times (*Quo ?A-0 ?A-*) 

(*Quo 1000000.0 ? L-By-D)) 

0.333)) 

We can get the value of Aspect Ratio from the user by 
issuing him a question . If he doesn't specify the value, 
a value of 8.0 is assumed. 

b. Area-zero is computed from the rule 

ST21 : (Area-zero ?A-0) -4- — — ( Overall-HTC ?U-0), 

(Log-Mean-Temp-Dif f ?LMTD), (Heat— rate ?Q ), 

(= ?A-0 (*Quo 2Q ( * times ? U-0 ? LMTD))) 

c. Overall-HTC has the following rule 

ST16 % (Overall-HTC ?U-0) 4 — ( Tube-side-fluid ?Tube-fluid) 

(Shell-side-fluid ? shell-fluid), 

(Fluid-type ?Tube-fluid ?Type-T), 

(Fluid-type ?Shell-fluid ?Type-S), 
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(Radius-Inside 9 R-l), (Radius-Outside ?R-0), 
(App-film-foul 9 Type-S 9 H-S 9 R-F-S), 

(App-film-foul 9 Type-t ?H-T 9R-F-T), 

(= ?U-0 

(*Quo 1.0 (+ (*Quo ?1.0 ?H-S) ? R-F-S 

(*times (*Quo 0.0025 ? wall-cond) 

(*Quo (*Times 2.0 ? R-0) 

(*Plus ? R-0 ?R-I ) ) ) 
(*times (* Plus 9R-F-T (*Quo 1.0 9H-T)) 
(*Quo 9 R-0 9R-I ) ) ) ) ) 

d. The shell side fluid is determined by 

ST3: (Shell-side-fluid 9Name- c ) < (fluid-name cold-fluid 

9 name- c ) , 

(Fluid-name hot fluid 9 name -h), 

(Tube-side-fluid 9 name-h) 

ST4 : (Shell-side-fluid 9name-h) ■ « — 

(Fluid-name Cold-fluid 9 name-c), 

(Fluid-name hot-fluid 9 name -h) ( Tube-side-fluid 
9 name-C). 

The shell-side-fluid has 2 rules, ST3 and ST4. 

The queries for fluid-name are satisfied by asking question. 

e. The tube-side-flui * has the following rules; 

ST5; (Tube-side-fluid 9 name-h) < — 

(Fluid-name hot-fluid 9name-h), (corrosive 

9name-h 9X), 

(iS-yes 9X). 
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ST7 : (Tube-side-fluid ?name-h) «e— 

(Fluid-name hot-fluid ?name-h), (fouling 

?name-h ?X) , 

(is-yes ?X) 

ST9: (Tube-side-fluid ?name-h) < — 

(Fluid-name hot-fluid ?name-h), 

(Fluid-name cold-fluid ?name-c) 

(Mass-flow-rate-cold ?name -c ?X), 

(Mass-flow-rate-hot ?name~ h ?Y) 

(> ?Y ?X) 

and similar rules for the cold-fluid. These rules are checked 
one by one. The rule ST5 states that if the hot fluid is 
corrosive, then it should be placed on the tube-side. The 
predicate is-yes checks whether the argument is 'yes' or not . 

Rule ST7 states if the hot fluid is fouling, then 
it should be placed on the tube side. The rule ST9 states 
that the tube-side has the fluid with the higher mass -flow- 
rate. Similar rules .exist for the cold-fluid. It can be 
seen that the either of the rule ST9 or STIO matches since some 
fluid will have mass-flow-rate higher than the other. At this 
step we infer the tube-side-fluid. When step c is complete, 
the new query is step d. This is also inferred since the 
tube-side-fluid is known. We then arrive at step c. The first 
two queries of the Antecedent are know. 3rd and 4th queries 
are satisfied by asking questions. This question asks the 
user to classify the hot and cold fluids in any one of the 
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given type. Queries 5 and 6 are satisfied by the rules, 
which require tube diameter and tube wall thickness. Tube 
Diameter and Tube-wall thickness are known from the question. 
Queries 7 and 8 are inferred from facts. The app— film— foul 
predicate has values of the approximate heat transfer coe- 
fficient and fouling factor for different types of fluids, 
stored as facts . After getting the various parameters, U Q is 
computed. The next query is Log-mean-temp-diff at state b. 

f. Log-men-temp-diff has the following rule. 

ST23: (Log-men-temp-diff ?LMTD) < — — 


(Fluid-name 

hot-fluid 

?name- 

-h). 


(Fluid-name 

cold-fluid ?name 

-c). 


(inlet-temp. 

?name-h 

?TH-I) , 



(Inlet-temp. 

?name-c 

?TC-I ) , 



(Outlet-temp. 

?name-h 

? TH-0 ) , 



(Outlet-temp. 

?name-c 

7TC-0), 



(= 7LMTD (*Quo (*Dif 

(*Dif 

7TH-I 

7TC-0 ) 



(*Dif 

7TH-0 

7TC-I ) ) 


(Log 

(*Quo 

(*Dif 

7TH-I 7TC-0) 




(*Dif 

7TH-0 7TC— I ) ) ) ) ) 


The fluid-name, Inlet-temp and outlet-temp are satisfied by 
asking questions. LMTD is computed from the equation. The 
next query at stage b is heat-rate • 



57 


g. The rule for heat-rate is 

ST22 : (heat-rate ?Q) ^ 

(Fluid-name cold-fluid ?name -c ) , 

(Mass-flow-rate-cold ?name-c ?M-c), 

(Specific-heat ?name-c ?Cp-c ) 

(Inlet-temp. ?name-c 7TC-I), 

(Outlet-temp. ?name-c 7TC-0), 

(= ?Q (* 7M-c ?CP-c (*Dif 7 TC~ I 7TC-0))) 

The first 5 queries are satisfied by asking questions to 
the user, in the previous steps. At this stage, the variables 
get the values directly from the ALIST. The heat-rate is 
computed from the equation given. Remembering that at stage 
b, we had to satisfy the subqueries at stage f and g, we 
observe that the Area-zero can be computed by the given 
equation, concluding step b. 


We move upwards at stage a. The first query is 
already satisfied at stage b. The next query is (Area-* 

7 A-*). 

h. Area-* has the following rule 


ST19 : (Area-* ?A-*) 

(Pitch-layout ? Type), 

(Constant-of-pitch-layout ?Type 7CPL), 
(Tube-diameter 7D-T), ( Tube-pitch 7L-TP), 

(= 7A-* (* 2.45 (*Quo 1.0 7 CPL) 

(*Quo ?D— T (Square 7L-TP) ) ) ) 



START ) 


58 


* 



Fig. 4.1 Flow chart for approximate design of STHE . 
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The pitch-layout is decided by asking question pertaining 
to pressure-drop and mechanical-cleaning. Staggered-30 type 
of layout gives a higher pressure drop and mechanical cleaning 
is difficult. Staggered-45 has ease in mechanical-cleaning. 
Inline-90 is preferred for low-pressure-drop. The constant 
of-pitch-layout has facts stored for each one of the 3 pitch- 
layout-types. Tube-Diameter is known by asking question. The 
tube pitch is computed as 1.25 times the tube-diameter by 
Rule no. ST20. 

ST20 : (Tube-pitch 7L-TP) 

(Tube-Diameter ?D-T), (= 7L-TP (* times 1.25 7D-T)). 

Area-* is computed from the equation. We move to the next 
query ( Aspect-Ratio 7 L-By— D) at the stage a. Aspect-Ratio 
is known by issuing a question. Hence the shell— diameter at 
stage a is inferred. 

j. Once the shell-diameter is inferred, the shell-length 
is computed by rule 

ST145 : (Shell-length 7L-TA) < ■ 

(Aspect-ratio ?L-By-D) , (Shell-Diameter 7 D- s ), 

(= 7 L-TA (*times ?L-By-D 7D-S)). 

The above steps can be easily summarised in the flow-chart 
f ig.4.1. 

Proceeding in the same way, the further rules 
were written. The remaining steps are described in short 
in the following sections. The detailed flow chart for the 
entire design process is given in Fig. 4.3. 



Property data 


(start) 

Jl 

Approximate 

design 
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Specification 
8. limitations 



Fig. 4.2 Flow chart for the design of STHE . 
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4.1.2 The second step is the evaluation of geometric 
parameters ( more detailed than earlier). Geometric 
parameters mean various flow areas? centriangle of baffle 
cut etc. These parameters are required for determining the 
correction factors for heat transfer and pressure drop. The 
basic dimensions of the tube are supplied by the user. The 
shell dimension are computed in the step 1 and at step 2, 
we have 3 further steps as 

(i) Determine shell and tube geometry. 

(ii) Determine baffle and bundle geometry. 

The evaluation of bundle geometry requires determining 
the bundle type. Bundle type is selected from the 
inlet condition of the fluids and as well as users 
requirement. 

Baffle cut is a function of shell diameter 
as well as central baffle spacing. The percentage 
of baffle cut determines the window correction factor. 

(iii) Determine the various flow areas for computing the 
correction factors. 

4.1.3 Determining various correction factors: 

There are methods available to compute the 
heat transfer coefficient and pressure drop for tube banks 
in pure cross flow. Since in the actual HE the flow is not 
pure cross flow? a measure of deviation from the cross flow 
is given by the correction factor. 
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4.1.4 Actual shell side heat tranfer coefficient can be 
easily computed once the correction factors and ideal tube 
bank heat transfer coefficient is known. 

4.1.5 Total shell side pressure drop is composed of 3 
distinct parts and is computed as a sum of the three pre- 
ssure drops (a) pressure drop in baffle windows (b) press- 
ure drop in cross flow and (c) pressure drop in end zones. 

4.1.6 Actual tube side pressure drop involves computing 
Reynolds number and Viscosity correction factor as the main 
steps. The tube side Reynolds number depends upon the 
velocity of the tube side fluid which in turn depends upon 
the no of tube side passes. 

4.1.7 Tube side pressure drop is computed after computing 
the Reynold number and the firction factor. 

4.1.8 Comparision of the results and iteration is the step 

which deals with modifications in the design. The process 
from step 1 to step 7 can be termed as one cycle. At this 
stage the user decides whether the design is acceptable or 
not. The design for which the heat transfer area 

required is less than the area available within the permi- 
ssible pressure drop is an acceptable design. The greater 
the difference between area required and area available, or 
between the pressure drop permissible and the actual pressure 
drop, poorer is the design. Such a design needs modification 




Fig. 4.3 Tree structure for the design of STHE . 
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Fig. a 


A - Value of the predicate argument (s) is known 
Q - Question 
F - Fact 
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4.2 GENERATING THE QUERY TREE s 

The query tree generated is shown in Fig. 4.2. 
Each block represents a sub query tree, details of which 
are shown in Fig. 4.4 . The blocks shown in Fig. 4.3 and 
4.4 are same except for the fact that Fig. 4.4 shows a 
detailed tree structure. 

It may be pointed out that the flow chart for 
the design will be the same but the tree structure depends 
totally on the structuring of the rules and the users 
responses. Hence it is a dynamic structure and may vary 
from session to session. In the present system, the queries 
are to be satisfied from the left to right. 

We will examine in detail a query tree for 
determining the bundle diameter. Bundle diameter is calcu- 
lated as the difference between shell diameter and the 
bundle-shell clearance. The bundle shell clearence is a 
function of shell-diameter and the bundle type. Bundle 
type depends upon the fluid properties like fouling, Inlet- 
pressure, non-lethal characteristics and mechanical cleaning 
frequency. Hence it can be seen that the query (Bundle- 
Diameter ?X) generates question like fouling characteristics, 
slightly greater diameter, non-lethal characteristics and 
slight leakage. 


In similar manner each of the block can be 
analysed to its detailed tree as shown in Fig. 4.3. The 
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symbol & denotes that the query for that particular 
predicate is satisfied at that stage or any stage previous 
to that. 

4.3 STRUCTURING OF THE RULES: 

As explained in Chapter 2, this system is a 
backward driven system. When a goal is given, the system 
tries to satisfy that goal, by scanning through the facts, 
rules, etc. 

A detailed rule heirarchy structure is shown 

in Fig. 4. 4. 

The top level predicate is 'complete'. To 
satisfy this goal, 3 subgoals should be satisfied viz. 
Area-req, Total-shell-side-spressure-drop and tube-side 
pressure drop. Area-req further requires 4 goals 
and so on. The order of goals to be satisfied is from 
left to right i.e. the left goal is to be satisfied first 
then the next one and so on. 

It is observed that some of the goals are 
required to be satisfied at many places. The context is 
declared to be true for all the predicates except computa- 
tional predicates. This means that all the inference is 
stored as a fact and when it is required next time in the 
subsequent rules, it is directly inferred from the fact 
itself thus saving the laborious task of applying the 
rule again and again. 
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This feature of storing the inference poses one 

problem. When a predicate is inferred, at some stage we may 
want the value of-’ this predicate arguments to be changed. 

Such a thing is possible by defining it as an action predicate 
or defining a predicate with another name. The process of res- 
eting the database is adopted here. 

The predicates whose value of the argument needs to 
be changed are identified and listed. The order of operat- 
ion for these predicates is changed as- 

1. Compute-pred 

2. Action-pred 

3. Facts 

4. Question 

5. Rules 

6. Question-proc . 

This differs from the regular order of operation as mentioned 
in Section 2 in the sense that the question is asked before 
the rule. 

So initially, when an approximate value is to be 
determined, the user has to say 'DONTKNfOW* so that the 
rule which computes the value is applied. 

After the preliminary round of rating is over, in the 
next round the properties are * Reset* to the orignal i.e. 
they lose all the newly generated facts and the system is 
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again ready to apply the same rules. This feature of asking 
the question is also useful when the user knows some piece of 
information and wishes to override the procedure that exists 
for computing that value. For e.g. if he already knows the 
tube-side heat transfer coefficient, he can give the value 
and override the usual procedure of computing the tube-side 
heat transfer coefficient 

4.4 DIFFERENT PREDICATES USED: 

The computational predicates take their argume- 
nts and return True (T) or false (Nil) according to the con- 
dition e.g. <,>, =, LE,GE etc. ( LE = less than or equal to 
and GE = Greater than or equal to.). Most of the predicates 
defined for this system take one argument which is a numerical 
value e.g. Heat-rate, no-of-tubes, shell-diameter. Such predi- 
cates essentially have a rule, and most of the time they are 
enf erred from a rule. They can get their value from the que- 
stion also as explained in the previous section. 

Other types of predicates are those which have 
more than one argument. The need for defining such predica- 
tes arises because of change in references and requirement 
of storing many properties. In the initial stage when the 
shell side and tube-side fluids are not decided, questions 
are asked by referring to the fluids as hot fluids and cold- 
fluids. All the inference till this stage is made by refe- 
rring the streams as hot and cold. Once the fluids on shell 



80 


and tube side are determined, they are referred by the 
predicates shell-side-fluid and tube-side fluid so (Inlet- 
temp. ?name ? T-l) is a predicate having name of the fluid 
as its first argument and the value as its second argument. 

The first argument can be had from either of the two predica- 
tes as explained - 

(fluid-name hot-fluid ? name) or 

(shell-side-fluid ? name) 

Another example is 

(App-f ilm-f oul ?type ? HTC ?Fouling-coeff ) 

This predicate stores the values of Heat-transfer coefficient 
and Fouling coefficients for different types of fluids. 

4.5 KINDS OF RULES: 

The rules used here are mostly of the type 
which use an expression for assigning value to a variable 
argument of the predicate . 

A rule of the form 

((( Area-req ? A-req) 

(( Overall-HTC-ac tual ? U-0) 

( LMTD - correction ? F) 

( Heat— Rate ? Q— 0) 

(= ?A-Req (* Quo ?Q-0 (* ? U-0 ?LMTD ?F)))) 

ST 133)) 

is equivalent to the expression 

Q o 

Areq = U q (LMTD)F 
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This can .be asserted as a logic expression as (Area-req 
?X) — ( overall-HTC-actual ? U-O), 

(Heat-rate ?Q-0), (LMTD-Correction ?F), 

( = ?X (*Quo ?Q-0 (* ?U-0 ? LMTD ? F))), 

The second type of rule is one which takes purely logical 
decisions. Consider the rule 

((( Shell-side-fluid ? name -C) 

(( Fluid-name cold-fluid ? name-c) 

(Fluid-name hot-fluid ?name -h) 
(Tube-side-fluid ? name - h) 

ST 3)) 

here this rule cannot be expressed as a mathematical state- 
ment, but it takes a logical decision that if it can be 
proved that there is hot-fluid on the Tube-side, then it can 
be inferred that the other (cold-fluid) is on the shell-side. 

4.6 ENCODING A NEW RULE; 

Here we will talk about creating new predicates 
and writing a rule. Consider the following statements 
' The decision of tube-side, shell-side or direct condensa- 
tion depends upon the following: 

1. High pressures are best inside tubes 

2. Shell side has a relatively low pressure drop. 

3. Corrosive vapours should be on the tube side. 

4. If the operating temperature is very high, use direct- 
condensation. 
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5. If the condensate can freze during the heat transfer 
process, shell side condensation is preferred. 

6. When condensing multicomponent mixtures having a sub- 
stantial boiling or dew-point range, or when there are 
soluble gases present, it is necessary to control the 
condensate and vapour flow so as to enable the low boi- 
lers to condense or when stripping to prevent their 
condensation or absorption, use tube side-condensation. 

7. Fouling vapours should be best placed on the tube 
side. 

8. Venting ;’If non-condensables are to be vented, use tube 
side condensation’ 


The first step is to arrange these statements in 
order of importance. This is the first and the most import- 
ant step towards writing the rules. 


Along with this, select the important parameters 
and call them as predicates. So we have- 


Order of Serial 

importance number 


Predicates 


1 

2 

3 


3 (corrosive <name> < yes/no>) 

5 (freezing <name> <yes/no> ) 

6 (multi-component-vapours <yes/no > ) 
or (soluble-gas-present <yes/no>) 
and (control-condensate-vapour flow 

<yes/no>) 
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Order of Serial Predicates 

importance number 


4 

5 

6 

7 

8 


1 (inlet-pressure ? name ? P-l) 

(< ?P-I 2000.0 ) 

4 (inlet-temp ? name ? I-I) 

(> ?T-I 300.00) 

8 (venting-noncondensables <yes/no>) 

2 (Pressure-drop-max-low <yes/no>) 

7 (fouling <name> < yes/no>) 


now using the syntax, rules can be written. One can have 
predicate of the type- 

condensation tubeside) or 

(Tube-side-condensation) The latter one is used 

here 

(Def prop Tube-side-condensation 
( ( (Tube-side-condensation) 

((Fluid-name hot fluid ?name -h ) 

(Corrosive ? name-h ?X ) 

( = ? X YES)) 

No 1 )) 

Rules ) 
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CHAPTER 5 

RESULTS AND DISCUSSION 

5.1 USING THE PACKAGE: 

The various steps are explained in sect. 3.1. 

The user is expected to identify the problem and select a 
STHE as a solution. This system does all the steps in the 
....... block shown in Fig. 5.1. While using this system, 

the user is expected to answer questions. The data are 
supplied essentially by question and answer process. 

Questions asked can be classified into 2 catego- 
ries . 

1. Mandatory questions and 

2. Optional questions 

1. Mandatory Questions: 

The user may have a requirement of heat exchange 
between any two fluids. Since the fluids undergoing heat 
exchange are best known to the user and so also are the 
properties, the user should supply them to the system when 
they are asked. For mandatory questions, if the user types 
’don't know' there is no way the system can get this value 
and hence the attempt fails. For example a typical mandatory 
question is 'give the inlet temperature of hot fluid in 
degrees Celsius ' • All the mandatory questions are self 
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USER 
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required 
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Fig.5.1 Methodology for using the system 















explanatory, and they demand certain values. There is no 
information stored for them ( See Table 5.1). 

2. Optional Questions z 

Some questions help as guide lines in selecting 
the bundle type, fluid side etc. If all the questions are 
not successful, there is still some decision taken, but it 
becomes a conservative or extra safe, decision. Answering 
the optional questions helps in getting the HE design pre- 
cisely to the user needs, and so user should attempt to ans- 
wer these questions. These questions are interdependent 
i.e. a question may be asked only if the user fails to answer 
the previous question : 

Thedatabase for the Rule based system for heat 
exchangers is stored in 2 files. 

FIXED AND ITER 

FIXED essentially consists of the parameters which will not 
be changed in the process of redesigning. They are essenti- 
ally the fluid properties and characteristics. ITER contains 
the bulk of the Rule base and is the one that rates the heat 
exchanger. Tho system consists of 148 Rules, 35 questions and 
explanations. The user is expected to fill up the tables I 
and II to facilate the use of this system. Table I contains 
the essential questions which he should know and table II 
contains some questions which may be required to be answered . 

The complete design of a heat exchanger is an 
iterative process and can be best handled by the user himself. 



A. After loading the files FIXED and ITER the user can 
start the system by typing Design. Before answering any 
question, he should read the instructions for answering the 
question carefully. If the user does not know the answer 
for the question, he should first get the explanation by 
typing 1 What * and see what are the alternatives he has got. 

At the end, a short comparative table is printed 
which gives the values required for comparison. 

This table is shown in the sample run . The user 
should examine the output and decide whether or not the design 
is acceptable. The following should be taken into consider- 
ation. 

1. Area required for heat transfer should in no case be 
greater than the area available for heat transfer. 

2. Area required may be smaller than the area available, 
but this descision can be made by the user himself. Larger 
difference is advantageous from reliability point of view, 
but is not economical. 

3. The required area and available area may be acce- 
ptable, but not the pressure drop on the shell and tube 
sides. The initial design may have the lowest pressure drop 
which will result in a very expensive design of the HE. 

The user should compare the actual pressure drop and 
the permissible pressure drop. A rough guide for determining 
the maximum permissible pressure drop is given below, (Lord, 
1970). 



Inlet pressure 
kPa 

Maximum permissible 
pressure drop 
kPa 

Sub atmospheric 

1/lOth of absolute pressure 

6.8 kPa to 68.9 kPa 
(gauge) 

1/2 of operating gauge 
pressure. 

68.9 kPa and higher 

34.5 kPa or higher 


A good design is the one which utilizes the maximum 
permissible pressure drop, and has the highest heat transfer 
coefficient. 


It may be noted at this point that the pressure drop, 
Reynolds number and Heat transfer area are mutually interdepent 
Heat transfer coefficient can be encreased by increasing the 
Reynolds number, but thi r increase in Reynolds number increases 
the pressure drop also. 

It is therefore, not possible to predict the exact 
behaviour of the heat transfer coefficient, pressure drop, 
and Area of heat transfer required simultaneously. A trial 
and error solution is the only way. 

4, Baffle leakage correction Factor for Heat Transfer (J^): 

A good design should have its value not less than 
0.6. It can be increased by increasing the cross flow area 
at shell centre or by decreasing the number of baffles. The 
initial design has the maximum value can attain and hence 



while iterating, case should be taken to see that it does 
not fall below 0.6. 
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5. Total Correction Factor; 

The product of the all 3 correction factors should 
be >_ 0.5 for a well designed H.E. and should never be < 0.4. 

B. Iteration; 

When the user is not satisfied, the system automati- 
cally iterates. The user should not down the values from the 
comparison table and take the decisions as explained in the 
previous^ section. He should write down what changes he wish- 
es to make under to corresponding columns ( see Table 5.3 ). 

The user can also save the previous session. 

During iteration only sleeted questions are asked 
and the user should answer them if he wishes to change the 
values. 

Many question are interdependent and the user should 
answer - those questions only e.g. The tube side passes and 
flow velocity intubes. Flow velocity inside the tubes in 
required to determine the tube side passes, but if the user can 
manipulate the number of tubeside passes himself, the 
question regarding tubeside velocity will not be asked. 

After this short session of questions, the same 


table of results will be displayed for the user to take his 
decision. 
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mandatory 


Question 


Name of the fluid 


Hot fluid Cold fluid 


Specific heat j/K g K 

Inlet temperature °C 

Outlet temperature °C 

Mass flow rate Kg/Sec 

(Liquids only) C P 

absolute viscosity at 
average temperature 

(Liquids only) C p 

absolute viscosity at 
some other temperature 


value value 
value value 

value value 

value 

value value 

value value 


Kg/m 3 


Density 


value 


value 
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TABLE 5.2 
ITERATIVE PROCEDURE 


SI. 

No. 

Parameter 

Units 

Itration no 

3 2 3 4 

1 . 

Heat transfer area required 

m 2 


2. 

Actual Heat Transfer area 
available 

m 2 


3. 

Shell side Re 



4. 

Tube side Re 



5. 

Shell side fluid 



6. 

Tube side fluid 



7. 

Film coefficient tube 
side 

W/m 2 K 


8. 

Film coefficient shell side 

W/m 2 K 


9. 

No of Baffles 



10. 

No of Tubes 



11. 

Overall Heat Transfer 

Coefficient 

W/m 2 K 


12. 

Total shell side pressure drop 

kPa 


13. 

Total Tube side pressure drop 

kPa 


i * 

0 

No of shell side pases 

kPa 
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It is observed that the solution is obtained faster 
if the actual heat transfer coefficient value from the previ- 
ous round is supplied, instead of supplying the area of heat 
transfer. 

After it is decided that there will not be any 
changes in the design, the final output is printed. 

Results are listed in the form of a table and are 
devided in 3 sections. They are presented at the end of 
Chapter 5. 

5.2 SAMPLE SESSIONS; 

Three sample sessions are presented here to illustrate 
how the expert system will respond if the required questions 
are not answered or answered in a typical manner ( see Appendix] 

Sample sessions I and II refer only to the mandatory 
question. If these questions are not answered, the system just 
fails. If these questions are answered in a different 

manner, it may result in an entirely different arrangement 
of the heat exchanger. For example if the user has been asked, 
"Is the hot gas significantly corrosive?' and the answer given 
is 'YES', this fluid will be placed on the tube side and the 
other fluid on the shell side. If the answer was *N0*, other 
factors will decide which fluid will be on the tube side and 
which one on the shell side. 

Sample session III gives an idea of a normal session 
to the reader. It starts with the system asking questions. The 
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user responds to the questions by seeking various explanations, 
etc. The session shows iteration for 3 times. In each case 
the users objective was to match the actual heat transfer area 
required and the area available for heat transfer. It can be 
noted that the user has to answer very few questions while 
iterating since all the other properties are stored as facts. 

The user has the facility of storing and recording the session. 

5.3 CONCLUSIONS: 

An 'Expert System' based upon the logic programming 
for the thermal design of single phase flow STHE has been deve- 
loped. The system is highly interactive which gathers informa- 
tion about the problem domain by asking the user relevant 
questions and giving him the necessary help and explanation (s ) 
if required. The system can be used pretty easily by a comm- 
ercial user for design purposes. It can also be used for 
learning heat exchanger design or teaching the same to a begi- 
nner. 

The detailed design procedure for STHE which consists 
of the shell and tube side designs, has been studied. Of all 
the (design) procedures available for the shell side analysis, 
modified Bell-Delaware method has been used. As the flow in 
the actual STHE is a non-ideal cross flow, the relevant corr- 
ection factors which can take this deviation into consideration 
have been calculated. Kern's method has been used for the 
tube side analysis. 
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The complete design procedure of STHE has been 
transformed into a knowledge base comprising of rules , 
facts, questions and explanations, which gives rise to 
a query tree indicating the flow of control. 

The expert system presented here can be modified to 
be able to add a new rule, change any of the value already 

given and obtain information about various parameters at the 

intermediate stage. It can be extended to incorporate the rules 

for the design of other types of heat exchangers. 
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